Multiscale Dynamics of Solar Magnetic Structures 
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ABSTRACT 

Multiscale topological complexity of solar magnetic field is among the primary factors control- 
ling energy release in the corona, including associated processes in the photospheric and chromo- 
spheric boundaries. We present a new approach for analyzing multiscale behavior of the photo- 
spheric magnetic flux underlying this dynamics as depicted by a sequence of high-resolution solar 
magnetograms. The approach involves two basic processing steps: ( 1) identification of tim ing 
and location of magnetic flux origin and demise events (as defined bv lDeForest et ail ( 2007|) ) by 
tracking spatiotcmporal evolution of unipolar and bipolar photospheric regions, and (2) analysis 
of collective behavior of the detected magnetic events using a generalized version of Grassberger 
- Procaccia correlation integral algorithm. The scale-free nature of the developed algorithms 
makes it possible to characterize the dynamics of the photospheric network across a wide range 
of distances and relaxation times. Three types of photospheric conditions are considered to test 
the method: a quiet photosphere, a solar active region (NOAA 1 0365) in a quiescent non-flaring 
state, and the same active region during a period of M-class flares. The results obtained show 
(1) the presence of a topologically complex asymmetrically fragmented magnetic network in the 
quiet photosphere driven by meso- and supergranulation, (2) the formation of non-potential mag- 
netic structures with complex polarity separation lines inside the active region, and (3) statistical 
signatures of canceling bipolar magnetic structures coinciding with flaring activity in the active 
region. Each of these effects can represent an unstable magnetic configuration acting as an energy 
source for coronal dissipation and heating. 

Subject headings: Sun: activity - Sun: magnetic fields - complexity 



1. Introduction 

The solar magnetic field is often seen as 
the underlying cause of solar activity and the 
primary agent coupling the Sun to the helio- 
sphere. High-resolution observations show that 
it possesses an incredible level of complexity 
across all relevant spatial a nd tempora l scale s 



(e.g.. ISchriiver et ai]([l992h: iBalke et all (fl993[> : 

hvnl 



Lawrence et al. Cad avid et al.l 

Abramenko et al.l (120021) :lAbramenko fc Yurchvshvr. 



(2010a); Georgoulis 



(2005): lMcAteer et al.l ([20101 ); 



Parnell et al.l (|2009l )). Understanding the origin 



and the physical significance of this complexity 
is among the most fundamental unresolved prob- 
lems of contemporary solar science. Its further 
analysis requires accurate quantification of cross- 
scale plasma processes involved in the formation 
and evolution of the observed multiscale magnetic 
patterns. 

Solar magnetic complexity is deeply rooted in 
the life cycle of coronal perturbations. Coro- 
nal magnetic field is known to accumulate pro- 
gressively large amounts of magnetic energy and 
helicity before producing explosive phenomena 
such as flares and coronal mass ejections (see 
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e.g. iSolanki et all (|2006l) ). The emergence o : 
complex domains of magnetic connectivity ac- 
companying this process is instrumental as it 
makes th e magnetic free e n ergy available for dis- 
sipation dAntiocho s 19981 : Schriiver et al. 2005; 



Edmondson et al. 2010() . The initial build-up 



of the magnetic complexity can be driven by 
the injection o f new magnetic structures (e.g., 
iFalconer et al. (|2008l): lAbramenko fc Yurchyshynj 
(|2010ah : lMcAteer et all (<2010l ): IConlon et al.1 hOlti ). 
or by the fragmentation of the existing flux sub- 
ject to turbulent footpoi nt shifting both in solar 



active regions (see e.g. Abramenko fc Longcopd 



(120051)') and in the qui et Sun (|Lamb et al.l (|2008): 



Chaouche et al.l (|2011 ) and refs therein). 

The resulting magnetic patterns featuring 
mixed-polarities regions with strongly sheared 
fields are prone to kinetic instabilities which can 
unfreeze the magnetic flux and start magnetic re- 
conncction. The reconnection brings the system 
to a lower energy state and reduces the mag- 
netic complexity l ocally but usually not glob- 
ally (jAschwandenl 120061 ) . Cooperative interac- 
tions between multiple reconnecting plasma tubes 
may lead to the second generation of multi- 
scale phenomena at the coronal level, includin g 
self-organized criticality ( Lu fc Hamilton! Il99l[) . 
storms of nan oflares in individual multistranded 
coronal loops ( Kli mchukll2006[ IViall fc Klimchukl 
20 111 iMorales fc Charbonneaul 120081 ). and cas- 
cades of in stabilities involving multiple recqnnect - 
ing loops (|Hughes et al.l 120031: iTorok et al.ll201ll) . 
The distortions in the topology of open field re- 
gions caused by some of these effects can be a 
major factor controlli ng both slowly varying and 



transient solar wind ( Antiochos et al.l 12007). In 



the quiet Sun, the wide ranges of spatial and tem- 
poral scales manifested by the brightness of the 
soft X-ray and EUV corona are likely to reflect the 
multiscale geometry and d ynamics of the under - 
lying photospheric network ( Falconer et al.lfl998h . 
The importance of solar magnetic complexity in 
triggering coronal dissipation events and the pos- 
sibility of their self-organization into multiscale 
structures is confirme d by many theoretic a l mod - 
els (see, for instance. ICharbonneau et al. |200ll): 



Aschwanden fc Parnelll (|2002h : lAschwandenl (|201ll )) 



Gaining a deeper insight into the nature of so- 
lar magnetic complexity and its impact on coro- 
nal instabilities requires innovative data analysis 



tools focused on cause and effect links between 
two or more interconnected multiscale processes. 
There is a serious need to explore their cross- 
correlation properties, in addition to the autocor- 
relation properties studied before. Furthermore, 
the causality of solar activity is an inherently spa- 
tiotcmporal phenomenon which unfolds both in 
spatial and temporal domains of analysis. Un- 
til now, the spatiotemporal aspect of solar mag- 
netic complexity has not been given proper atten- 
tion in observational studies. Thus, for example, 
the algorithms of solar magnetic tracking com- 
monly used for analy zing the emergence, cancel- 
lation (jMartinl Il984j ). fragmentation and coales- 
cence dynamics of the photospheric network are 
essentially monoscale. The outputs of these al- 
gorithms and the interpretation of the obtained 
results depend, explicitly or implicitly, on ad hoc 
choi ces of spatial and tem poral correlation scales 
(see DeForest et al.l ( 20071 ) for a detailed review). 
To our knowledge, none of the methods proposed 
so far can simultaneously address short- and long- 
range properties of solar magnetic network within 
the same data analysis framework while treating 
both spatial and temporal scales as free param- 
eters. The available post-processing techniques, 
such as Fourier spectrograms, spatial correlation 
functions, or multifractal statistics, share the same 
methodological "blind spot" - they cannot ade- 
quately describe transient cross-scale coupling as- 
sociated with time-dependent nature of studied 
multiscale patterns. 

In this paper, we propose a new data analy- 
sis paradigm focused on multiscale spatiotempo- 
ral properties of the photospheric flux dynamics. 
As opposed to previously proposed methods for 
identifying correlated magnetic events, the devel- 
oped algorithms contain no a priori assumptions 
regarding the intrinsic spatial and temporal inter- 
action scales. By accepting this strategy, we take 
into account the fact that magnetic concentration 
regions, acting as nodes of the photospheric net- 
work, can communicate over arbitrary long dis- 
tances beyond the local neighborhood. We also 
assume that these communications can encompass 
a wide range of response times. 

Our technique involves two main processing 
steps. First, we determine timings and locations of 



origin and demise magnetic events ((DeForest et al 



20071 ) based of a direct spatiotemporal tracking of 
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unipolar and bipolar photospheric regions above a 
specified magnetic flux threshold. Second, we an- 
alyze pairwise correlations between the detected 
magnetic events using a generalized Grassberger - 
Procaccia correlation integral technique. By de- 
sign of the algorithms, any detected characteristic 
scale is an attribute of the studied data and not 
that of the method. 

This proposed data analysis framework is suc- 
cessfully tested on a sequence of high cadence 
SOHO MDI magnetograms representing a so- 
lar active region (NOAA 10365) and a quiet 
Sun. The application of spatiotemporal method- 
ology extends earlier statistical analyses of mul- 
tiscale geometr y of th e photospheric fi eld (e.i 



Lawrence etaH (|l993l ); iMcAteer etahl (|201(j )) 



and allows us to evaluate scaling characteristics of 
the dynamics of the photospheric flux not reported 
before. The results obtained reveal significant dif- 
ference in the network dynamics underlying flar- 
ing and non-flaring regimes of NOAA 10365, and 
deliver a detailed picture of spatiotemporal clus- 
tering in the adjacent quiet photosphere. 

The paper has the following structure. Section 
2 provides the description of the studied data, so- 
lar conditions, and applied statistical tools. Sec- 
tion 3 presents main results of our analysis and 
their qualitative interpretation in terms key mul- 
tiscale plasma processes operating in solar pho- 
tosphere and corona. Section 4 summarizes our 
findings and highlights some remaining problems. 

2. Data and methods 

2.1. Observational data and solar condi- 
tions 

We analyzed a sequence o f 816 high-resolution 



SOHO MDI magnetograms (jScherrer et al.l 11995) 



collected between 14:00 05/26/2003 and 04:35 
05/27/2003, average cadence time 64 seconds. 
The images were calibrated using standard Solar- 
Soft procedures and rebinned down to the 0.9x0.9 
Mm resolution to reduce the pixel noise, with the 
final image size 512x512 pixels. The studied field 
of view contains two adjacent spatial domains rep- 
resenting different levels of solar activity (Fig. [TJ). 

The quiet region to the north of the solar equa- 
tor (labeled by letter Q in Fig. [TJ exhibits a char- 
acteristic pattern of strongly fragmented magnetic 
flux loosely organized into an irregular cellular 



structure typical of a quiet Sun. As we discuss be- 
low, this structure is associated with photospheric 
meso- and supergranulation, and it demonstrates 
a significantly geometric asymmetry of positive 
and negative magnetic polarities precluding di- 
rect identification of canceling bipolar pairs. The 
region operates in a statistically steady state as 
evidenced by the similarity of scaling properties 
of origin and demise dynamics across the studied 
range of scales. 

The southern portion of the studied field of view 
(labeled with A in the figure) is dominated by the 
NOAA 10365 active region which first appeared 
on the visible solar disk on 19 May 2003 at he- 
liographic latitude ~ —5°. A few days later, a 
new bipolar structure formed inside the active re- 
gion. By merging with the decayed preexisting 
flux, the new dipol e developed a comple x polar- 
ity separation line ([Verbeeck et al l l201lh . Start- 
ing from 25 May, the total magnetic flux of NOAA 
10365 was increasing steadily before reaching a 
plateau on 29 May 2003. The region gradually ro- 
tated counter clockwise with th e leading negative 
flux polarity (|Tian et al.1 1201 lb . The continuous 
injection of a new flux in the northern portion of 
NOAA 10356 created a 5— spot configuration with 
a positive region progr essively surrounded b y the 
emergent negative flux ( Chandra et alJ 20091 ) . The 
magnetic field pattern of the emerging bipole fea- 
tured two large elongated areas of opposite po- 
larity s uggesting a twisted flux tube of p ositive 
helicity (jTian et al.ll201ll IChae et alj|2004f ). 

The activity of the NOAA 10365 region under- 
went significant changes during the studied time 
interval. Between 15:00 05/26/2003 and 20:48 
05/26/2003, the active region generated multiple 
M-class flares, with the GOES 12 1.0-8.0 A flux re- 
maining above the 10~ 10 watts m -2 level. During 
20:48 05/26/2003 - 02:42 05/27/2003, the GOES 
flux stayed consistently below this level indicat- 
ing no noticeable flare activity. The two intervals 
were analyzed separately in order to compare mag- 
netic clustering effects in flaring and non-flaring 
regimes. 

2.2. Definition of magnetic events 

To determine timings and positions of magnetic 
events in the studied set of magnetograms, we 
i nvoked direct spat i otemp oral tracking approach 
(jUritskv et al ]l2007ll2010ah in which magnetic fea- 
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B<0 B>0 




Fig. 1 — SOHO MDI magnetogram (14:00 05/26/2003) showing the studied field of view (approx. 445x445 
Mm, centered at the solar disk center) including the subregions with quiet (active) conditions above (below) 
the dashed horizontal line representing solar equator. The NOAA 10365 active region is outlined by the 
dotted circular line. 
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tures living for more than one image sampling in- 
terval were treated as subvolumes of the three- 
dimensional space-time. The history of each fea- 
ture is represented by a connected set of image 
pixels containing an unsigned magnetic field above 
a specified detection threshold B t h- 

The idea of the threshold-based spatiotempo- 
ral tracking has been impl emented in the c ontex t 
of solar mag n etogr ams by iBerger fc Title] (|l996h : 
Berger et al.l (Il998h . The object tracking method 



ual image followed by cross-frame feature identifi 



cation within an image set, see e.g. iDeForest et al 



used in their study was based on measuring the 
centroid locations of the most compact and in- 
tense magnetic features. Here we invoke a some- 
what different approach which does not rely on the 
centroid location as a guiding center for feature 
tracking. Time-evolving magnetic concentrations 
are treated as true 3-dimensional objects, which 
makes it possible to study a more general class 
of magnetic structures whose positions cannot be 
adequately described by a centroid. 

In order to overcome the inherent memory lim- 
itations of standard cluster analysis algorithms 
such as, e.g., the Label Region function of the In- 
teractive Data Language (IDL), we have developed 
a new numerical scheme enabling fast decomposi- 
tion of a multidimensional data array into a set 
of connected pixels (clusters) while dramatically 
reducing the amount of stored information. 

The first step of our technique consists of build- 
ing a table of so-called "activations" - contiguous 
time intervals during which the condition \B\ > 
B t h is fulfilled in each of the image pixels. The 
second step is to find and label spatially con- 
nected clusters of the activations. This is done us- 
ing the breadth-first search principle which avoids 
backtracking of search trees representing individ- 
ual clusters. We find that it is important to con- 
sider all 8 nearest neighbors of each spatiotempo- 
ral node, including the diagonal neighbors, when 
identifying connected activations. Finally, the ta- 
ble of activations is rearranged according to the 
assigned cluster labels, in order to provide faster 
access to the detected structures. The output data 
array contains the complete information on the 
location, shape, and time evolution of all the de- 
tected magnetic concentration regions in the stud- 
ied 3D volume. 

It can be seen that the proposed algorithmic so- 
lution reverses the usual order of feature tracking 
steps (feature segmentation within in an individ- 



(|20071 )). Instead, we start with detecting temporal 
traces of features in individual photospheric loca- 
tions, and then consider their spatial adjacency. 
We found this algorithmic solution to be remark- 
ably fast and memory-efficient, with the size of 
the output data array being typically by two or- 
der of magnitudes smaller compared to the origi- 
nal data. The proposed scheme is obviously easy 
to parallelize. A more complete description of the 
developed algorithm and its numerical validation 
on high-Reynolds number magnetohydrodynami c 
turbulence can be found in (jUritskv et al.ll2010bl ) . 

Using the described tracking technique, we 
identified spatial and temporal coordinates of 
origin and demise magnetic events described by 
positive and negative line-of-sight magnetic polar- 
ity. The origin events (OE) represent newly de- 
tected magnetic concentration regions which just 
exceeded the threshold. The demise events (DE) 
occur when a previously found region becomes un- 
detectable, i.e when its flux content drops below 
B t h- This classification is analogous to the subdi- 
vis ion into " origin" and " demise" events described 
by IDeForest et al. ( 2007 ) although these authors 
used a different feature tracking method. It should 
be emphasized that our definition of events is not 
limited to the injection of new magnetic bipoles 
into the photosphere and their subsequent sub- 
mergence into the convection zone. In essence, we 
deal with an apparent dynamics associated with 
crossings of the detection threshold B t h by local 
magnetic field. While the resulting sets of events 
may contain a fraction of truly emerging and sub- 
merging magnetic bipoles, they also represents 
magnetic flux coalescence, fragmentation, cancel- 
lation, and other effects leading to spatiotemporal 
variability of the photospheric field. 

The photospheric velocities of positive an neg- 
ative magnetic spots inside the NOAA 10365 ac- 
tive region were in t he range 0.07 to 0.24 km/s 
(jChandra et al.ll2009l) . or less than 0.015 Mm per 
minute. For the time scales r < 20 min stud- 
ied here, this motion would produce a spatial dis- 
placement of < 0.2 Mm which is below the spatial 
resolution of SOHO MDI data. Based on this esti- 
mate, the magnetic flux configuration of the active 
region can be considered quasi-static in the photo- 
spheric plane, with most of the observed changes 
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owing to enhancement and reduction of local pho- 
tospheric magnetic flux rather than its horizontal 
convective motion. The same conclusion applies 
to the quiet region. Typically, the quiet Sun hor- 
i zontal speeds lie in the range 0.45 - 0.50 km/s 
(jShine et alj|2000l ); over a 20-minute interval, this 
motion is likely to be limited to a single MDI pixel 
and remain unresolved. 



To reduce the level of granulation noise ([DeForest et al 



20071 ). we focused on the magnetic features which 
lasted more than 5 minutes and encompassed a 
sufficiently large spatial area (> 4 Mm 2 ). The 
valid features were also required to have small 
initial and final linear sizes (< 3 Mm) ensuring 
accurate determination of their locations. The 
magnetic features which did not meet these qual- 
ity constraints as well as those truncated by the 
boundaries of the chosen field of view or observing 
interval were excluded from further analysis. 

2.3. The generalized correlation integral 
(GCI) analysis 

The structure of multiscale data sets such as 
the ones generated by solar magnetograms can be 
accurately described using th e family of order- q 
generalized fractal dimensions ( Grassbcrge r 1983): 



D q = 



d log e 
q — 1 d log e 



,9 = 1 



(1) 



otherwise. 



Here, n(e) is the number of non-empty hyper- 
cubes of a partitioning with length scale e, and 
Pi is the probability for a da ta point to fall into a 
cube i of the partitioning; see lSchuster fc Wolfram 
(|2005l ) for a pedagogical introduction to the con- 
cept. Among the dimensions defined by eq.(fT]), the 
Hausdorff fractal dimension (q — 0), the informa- 
tion dimension {q = 1), and the correlation dimen- 
sion (q = 2) are widely used. We base our analysis 
on the latter because it represents a second-order 
statistical moment which can be naturally con- 
verted into a cross-correlation measure of bipolar 
magnetic structures as shown below. 

For a set of point objects embedded in a multi- 
dimensional metric space, the correlation dimen- 
sion is commonly evaluated usi ng the Grassberger 

- Pro caccia correlation integral (jGrassberger fe Procaccial 
19831 ). The correlation integral (CI) is computed 



by counting the number C of pairs of set elements 
separated by the Euclidean distance < I. The 
shape of the C(£) dependence characterizes clus- 
tering effects in the studied set. For a random set 
with statistically independent element locations 
C(£) oc l d , where d is the embedding dimensions. 
For example, if the elements are uniformly scat- 
tered on a smooth surface (d = 2), the probability 
to find a neighbor within a given distance £ scales 
as the area of a circle with this radius, and hence 
C oc £ 2 . If the elements are clustered there exists 
an interval of scales in which C(£) oc £ a<d . The 
power-law exponent a of the correlation integral is 
an empirical proxy for the correlation dimension. 
If a is constant across a substantial range of scales 
(one ore more orders or magnitude), the studied 
set demonstrates self-similar or self-affine density 
variations involving voids and aggregates of vari- 
ous sizes. The discrepancy between the correlation 
dimension a and the embedding dimension d is a 
sensitive measure of this clustering behavior. 

The CI technique can be generalized to the case 
when the studied set contains a mixture of ele- 



ments of two distinct types (jKantzl 119941 ) . In our 
study, the two types of events are associated with 
positive and negative line-of-sight polarity mag- 
netic events (origin or demise) embedded in space 
and in time. The events are described by a com- 
bination of spatial coordinates r\ and timings ij, 
i = 1, .., N , where N = N + + N- is the total num- 
ber of positive and negative events. The subsets 
associated with OE and DE dynamics are treated 
independently. 

As the preparatory step, we compute the ma- 
trix L = {Lij} of pairwise distances L+j = | — r_j | 
and the vector T = {Ty} of pairwise time inter- 
vals Tij = \ti — tj\ between all the events, where 
i,j 6 We also define the binary vector 

P = {Pi} taking the value 1 (0) if the ith event 
is of positive (negative) type. Next, we construct 
the inclusion operator 



Ly)0(r-Ty) 



(2) 



in which 9 is the Heaviside step function. Note 
that ^ij(£,r) equals 1 if the two events i an j 
lie within spatial distance £ and their timing dif- 
ference does not exceed r, and takes the value 
otherwise. Finally, we define the following set of 
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generalized correlation integrals (GCIs): 



c_(£,r)=-l^x;(i-p i )(i-p i )*y^T) (4) 

C±(£, r)= (JV+) 1 (iV _ ) E ^(l-^)^ij(^ r) (5) 



Each of the GCI functions (C+, C_, and C±) 
approximates the probability of finding a pair of 
events of a particular type (respectively positive, 
negative, and bipolar) within the spatiotemporal 
neighborhood V = r oc £ d r in which is the 
volume of the ci-dimensional hypersphere (Fig. [2|. 
By construction, the autocorrelation integrals C+ 
and C- characterize the internal structure of pos- 
itive and negative subse ts of events and are anal- 
ogou s to the usual CI (jGrassberger fc Procaccia 
1983), whereas the cross-correlation integral C± 
describes interactions between positive and nega- 
tive events in bipolar magnetic pairs. The math- 
ematical definition of C± is similar to t he cross- 
correlation sum proposed bv lKantz (1994) for eval- 
uating an overlap of two chaotic attractors embed- 
ded in a reconstructed phase space of a nonlinear 
dynamical system. The primary difference is that 
here we deal with an essentially anisotropic met- 
ric space spanned by spatial and temporal coordi- 
nates of solar magnetic events. 

It is easy to show that randomly distributed 
events with statistically independent positions 
and occurrence times yield C oc V, VC £ 
{C + ,C-,C±}. Any departure from this "'null'" 
scaling regime indicates spatial and/or temporal 
correlations between the events. By adopting the 
scaling ansatz developed in th e theory of nonequi- 
libriu m critical phenomena ()Barabasi fc Stanley 
19951 ). we assume that C(£) and C(r) dependences 
can be factorized after a proper renormalization 
resulting in a set of spatial and temporal power- 
laws: 



C+ <x£ a +T p + 
C- oi£ a -T f} - 
C± (x£ a± T fj± 



(6) 
(7) 
(8) 
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Fig. 2. — Schematic illustration of the general- 
ized correlation integral (GCI) technique for d = 
2. Black (white) circles show positive (negative) 
events embedded in the d + 1-dimensional space- 
time. The cylinder represents a spatiotemporal 
neighborhood centered around the event outlined 
by the square and defined by a particular choice of 
£ and r. The outlined (positive) event has N + — 7 
unipolar and N± = 6 bipolar neighbors. The 
GCI calculation consists of an ensemble-averaging 
over all such neighborhoods centered successively 
at each of the events in the studied sets, for all 
relevant combinations of £ and r. 



The spatial autocorrelation dimensions a+ and a_ 
characterize geometric clustering of respectively 
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positive and negative events; the events of either 
type are clustered if the corresponding dimension 
is below d. Temporal autocorrelation dimensions 
f3 + and /3_ have the same meaning in the time do- 
main, with (3 < 1 reflecting clustering effects along 
the time axis (since the embedding dimension of 
time is f). 

As opposed to a+, a_, /3+ and /?_ dimensions 
describing unipolar correlations within the sets of 
positive and negative magnetic events, the cross- 
correlation dimensions a± and f3± arc intended to 
quantify mutual coupling between the two types of 
events. If negative and positive events are statisti- 
cally independent in space and in time, C± scales 
as oc r d r similarly to the autocorrelation integrals 
C+ and C_, and so a± = d and /3± = 1. For cor- 
related sets of bipolar events, the cross-correlation 
dimensions take values below or above the embed- 
ding dimension depending on the sign of the bipo- 
lar correlation (positive or negative, correspond- 
ingly). 

Asymptotically, the cross-correlation dimen- 
sions are independent of the in trinsic clust ering 
structure of the compared sets (|Kantdll994l ). In 
other words, if each of the unipolar sets is clustered 
but the sets are unrelated to each other, we expect 
that a± « d and j3± « 1 for sufficiently large N 
and negligible finite-size effects (so that t and r 
are much smaller than the largest relevant spatial 
and temporal scales). Fig. [3] illustrates this prop- 
erty of the GCI statistics on static sets of bipolar 
events constru cted by using t he stochastic Cantor 
set algorithm (|Turcottejll989T ) with the theoretical 
fractal dimension dF =log(3)/log(2). In panel (a), 
positive and negative polarity sets are generated 
by randomly sampling the same underlying Can- 
tor set. They have a statistically similar shape, 
and all three correlation integrals (C+, C_ and 
C±) collapse onto a single curve, with the log- log 
slope approaching dp as expected. In panel (b), 
the compared sets are created by sampling two 
independent runs of the Cantor algorithm. These 
sets have the same autocorrelation properties as 
in Fig|3][a) but must be uncorrelated with each 
other. Indeed, the measured a+ and a_ dimen- 
sions are close to the theoretical value dF, whereas 
a± approaches the embedding dimension d = 2 
confirming mutual statistical independence of the 
two sets. The GCI method leads to this correct 
conclusion despite the complex structure of the 
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Fig. 3. — Identifying auto- and cross-correlations 
in simulated clustered sets using the GCI ap- 
proach. The studied sets (blue and green dots) 
are generated by the same (a) and different (b) 
realizations of the underlying stochastic Cantor 
set (shown with gray color) with fractal dimen- 
sion log(3)/log(2) «1.158. Bottom panels (c) and 
(d) show auto- and cross-correlation integrals (eq. 
[3]- [5} computed for the pairs of sets in (a) and (b), 
correspondingly. The autocorrelation integrals of 
each set are shown with the matching color (blue 
or green). The cross-correlation integrals are in 
red. The code of the stochastic Cantor set is cour- 
tesy of G. Uritsky. 

compared sets generating false visual similarities 
at various scales. 

For the sake of brevity, the experimental part of 
this paper focuses on spatial scaling properties of 
magnetic events as reflected by a exponents com- 
puted at different r scales. The delay time is var- 
ied but treated as a parameter rather than as an 
independent scaling variable, with the main em- 
phasis placed on the multiscale geometry of unipo- 
lar and bipolar magnetic structures as well as the 
mechanisms of their cross-scale interaction and or- 
ganization. The analysis of j3 exponents and the 
associated temporal clustering is left for future re- 
search. 

To estimate a dimensions, we evaluated log-log 
slopes of the corresponding GCI functions using a 
least-square linear regression model built for sev- 
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eral selected ranges of spatial scales. In addition 
to this approach yielding a discrete set of a val- 
ues, we used the continuous scalogram technique 
introduced bv lUritskv et al. ( 2011 ) which enables 
classification of clustering regimes across the en- 
tire studied range of scales. The GCI scalograms 
are defined as 



a(£,r) 



d log C(1,t) 
d log I ' 



(9) 



where a G {a + ,a-,a±\ and C G {C+,C_,C±}. 
As in the case with the discrete GCI dimen- 
sions, the scalograms © were computed indepen- 
dently for the OE and DE populations of magnetic 
events. 

The statistics reported in the next section are 
for a single detection threshold B th = 150 G. Most 
of our conclusions holds for a wide range of B t h 
within 1-4 standard deviations above the mean 
magnetic field magnitude, with the exception of 
the cross-correlation dimension of the active re- 
gion exhibiting a more subtle Bth dependence as 
discussed in Section 3.3. 

3. Results and discussion 

3.1. Overview of scaling properties of the 
magnetic events 

Our analysis has shown that the cooperative 
behavior of the same and mixed polarity magnetic 
events depends strongly on their spatial and tem- 
poral separation as well as on the local solar ac- 
tivity level. It has been also found that bipolar 
events of either type (origin or demise) exhibit 
considerably different spatial correlation patterns 
compared to the corresponding type of the unipo- 
lar events. The difference is statistically signifi- 
cant and has a tendency to increase with coronal 
activity. 

Fig. [4]presents the combined results of the CGI 
analysis of OE and DE activity detected in the 
quiet and active photospheric regions. The left 
panels show the CGI functions (eq. - (O) eval- 
uated for two fixed time scales, r = min (dot- 
ted lines) and r = 10 min (solid lines). It can 
be seen that the shape of the autocorrelation in- 
tegrals characterizing unipolar magnetic clusters 
changes dramatically with the delay time. For 
r = lOmin, both C+ and C_ exhibit a distinct 



concave crossover (d 2 C/d£ 2 > 0) in the vicinity of 
t = £q ~ 15 Mm. As we discuss below, this spatial 
scale can be associated with the lower size limit of 
photospheric supergranules. 

The right panels in the same figure show the de- 
pendence of the GCI dimensions measured below 
(dashed lines) and above (solid lines) the crossover 
scale £q on the time scale r. The systematic de- 
crease of the autocorrelation dimensions a+ and 
a_ is likely a result of the P-mode contamination 
due to the five-minute Doppler oscillations leak- 
ing into the Zeeman signal and/or the granulation 
noise described by ro ughly the same time scale 
(jDeForest et al.l 120071 ). Each of these mechanisms 
can modulate the magnetic flux causing the event 
detection algorithm to loose tracks of weak mag- 
netic features with a field magnitude close B t h, 
and to rediscover them after one disturbance cy- 
cle. In either case, the increase of r beyond the 
5- minute time scale leads to better statistical rep- 
resentation of long-living photospheric structures. 
It is worth noticing that the cross-correlation in- 
tegrals C± (shown in red) are generally steeper 
than C-f and C_ , and unlike the autocorrelation 
integrals reveal no systematic dependence on r. 

Table [1] summarizes our measurements of the 
CGI dimensions for the quiet Sun as well as for 
the active region operating in the quiescent and 
flaring states. Sections 3.2 and 3.3 below provide 
a detailed discussion of the scaling properties char- 
acterizing these three regimes. 

3.2. Quiet Sun 

As can be seen from Fig. [4ja,b) an d Table [TJ 
unipolar dimensions a+ and a_ of the quiet pho- 
tosphere are systematically below the value 2, re- 
vealing clustering effects at the levels of meso- and 
supergranulation. At distances I < £o, these di- 
mensions are in the range to 1 for both origin 
and demise events. At larger spatial scales, the 
unipolar dimensions take values between 1 and 
2. The observed scaling crossover uncovers two 
distinct clustering regimes characterizing the ge- 
ometry of the photospheric flux below and above 
the supergranular scale £q. The small-scale dimen- 
sions are most likely to represent the tendency of 
the magnetic concentration regions to align along 
the boundaries of the photospheric supergranules. 
The non-integer a values indicate the presence of 
gaps and clusters of various sizes analogous to frac- 
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Fig. 4. — Results of the GCI analysis of origin and demise activity in quiet and active photospheric regions. 
Left panels: C+, C_ and C± plots for two time scales r = (10) min shown with dotted (solid) lines. 
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Table 1: Estimated values of GCI dimensions of origin (OE) and demise (DE) magnetic events in the quiet 
(Q) and active (A) regions. 

Region a + a_ a± a + a_ a± 



OE, t = min 

Q 1.66 ± 0.20 1.36 ± 0.06 2.65 ± 0.41 1.64 ± 0.01 1.19 ± 0.02 1.86 ± 0.01 

A 2.03 ± 0.27 2.01 ± 0.13 2.06 ± 0.18 1.46 ± 0.01 1.56 ± 0.01 2.24 ± 0.02 

DE, t = min 

Q 1.64 ± 0.13 2.01 ± 0.18 3.64 ± 0.39 1.64 ± 0.01 1.24 ± 0.02 2.47 ± 0.01 

A 1.30 ± 0.07 1.77 ± 0.05 1.25 ± 0.06 1.33 ± 0.01 1.51 ± 0.01 1.95 ± 0.02 

OE, t = 10 min 

Q 0.43 ± 0.02 0.72 ± 0.02 3.28 ± 0.13 1.40 ± 0.01 1.01 ± 0.01 2.25 ± 0.01 

A 0.78 ± 0.02 0.55 ± 0.02 2.02 ± 0.06 1.23 ± 0.01 1.43 ± 0.00 2.02 ± 0.01 

DE, t = 10 min 

Q 0.43 ± 0.01 0.63 ± 0.02 3.53 ± 0.15 1.46 ± 0.01 0.97 ± 0.00 2.32 ± 0.01 

A 0.86 ± 0.01 0.68 ± 0.02 1.84 ± 0.07 1.21 ± 0.01 1.41 ± 0.00 2.03 ± 0.01 

OE, t = 20 min 

Q 0.35 ± 0.02 0.59 ± 0.02 3.10 ± 0.11 1.38 ± 0.01 0.96 ± 0.01 2.26 ± 0.01 

A 0.72 ± 0.01 0.54 ± 0.02 1.94 ± 0.06 1.21 ± 0.01 1.40 ± 0.01 2.00 ± 0.00 

DE, t = 20 min 

Q 0.39 ± 0.01 0.57 ± 0.02 3.31 ± 0.12 1.38 ± 0.01 0.95 ± 0.01 2.33 ± 0.01 

A 0.72 ± 0.01 0.56 ± 0.02 1.70 ± 0.06 1.22 ± 0.01 1.37 ± 0.00 2.04 ± 0.01 
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tal clus tering in a on e-dimensional Cantor dust 
model (|Turcottdll989h . The large-scale cluster- 
ing at t > £o described by a_,a+ £ [1,2] re- 
flects a more isotropic spatial arrangement of solar 
magnetic flux. At these scales, the flux should be 
organized into quasi-two-dimensional supergranu- 
lar structures also containing multiscale voids and 
clusters. 

It is of interest to relate the obtained statistics 
with commonly discussed theoretical scenarios of 
photospheric supergranulation. The supergran- 
ular flow is known to co ntain a dominant hor - 
izontal component (e. g., iLeighton et al. ( 1962 ): 



Hathaway et al 



<|2000h ). Its leading role in the 
formation of the magnetic network outlining su- 
pergranula r cell boundarie s has b een intensively 
discu s sed (lHagenaar et ail Il997 ; Schriiver et al. 
11997k ISrikanth et al.ll2000() . In deterministic ex- 
planatory models, photospheric network appears 
as a result of forced rearrangement of small-scale 
surfa ce magnetic fi e lds by a large-scale horizontal 



flow (jMartini Il988t iLisle et allfeoOOt IWang et al 



1996). In stochastic models, the supergranulation 



network is generated by kinematic motion of ele- 
mentary flux tubes transported passively by the 
supergranular flow superposed with random dis 



placements produced by granulation ([Simon et al 



1995L 120011 ). by mutual advection of thermal down- 
flows taking place at random locations (jRastl 



2003), or as a self-organization of many interact- 



ing mag netic elements follo wing a random walk 
pattern (jCrouch et alj|2007t ). 

The deterministic models operate with a well- 
ordered supergranular flow which would result in 
integer a values both at small and large spatial 
scales, contrary to our observations. The second 
(stochastic) group of models seems to be in a qual- 
itative agreement with the correlation analysis re- 
sults. Our measurements show that the network 
is significantly fragmented at the supergranulation 
scale and exhibits an approximate co-alignment of 
magnetic flux elements with the lines of merging 
of the adjacent supergranular flow cells. The frac- 
tional values of a_ and a+ exponents at both 
small and large scales can be a footprint of a 
stochastic generating process, possibly related to 
the aggregation, fragmentation, and cancellation 
of randomly shuffled magnetic concentra tion re- 
gions as proposed bv lCrouch et al. ( 2007 ). 

To some extent, our estimates can be affected 



by the Doppler noise due to the flow motion at the 
supergranulation scale. In contrast to the short- 
scale noise sources mentioned in section 3.1, the 
supergranulation flows are capable of injecting ad- 
ditive signals into MDI magnetograms over a wide 
range of scales (see e.g. Meunier et al. ( 2010t ) and 
refs therein). However, we expect this mechanism 
to have a limited effect on our results. The rms 
values of the erroneous magnetic field caused by 
Doppler noise are known to increase with mean 
magnetic field. The error should be quite small 
for the magnetic field magnitudes comparable with 
the chosen d etection threshold B tn- The simula- 
tion study bv lSettele et~aTl (|2002l ) predicts SOHO 
MDI rms errors to be less than 10 G for a relatively 
strong mean signal of 800 G subject to 5-minutc 
velocity oscillations with amplitude up to ~ 300 
m/s. Extrapolating their estimates to B t h = 150 
G yields an uncertainty of ~ 1 G which can be 
safely neglected. This conclusion applies to both 
quiet Sun and AR 10365 active region since we ap- 
plied the same detection threshold on both sides 
of the solar equator. 

Furthermore, the supergranulation Doppler 
noise exhibits strong dependence on the helio- 
centric angle from disk center. At large heliocen- 
tric angles, the line-of-sight motion is dominated 
by the horizontal component ( 1 00-300 m/s, see 
e.g. estimates by IWelsch et al. ( 2012 ) based on 
local correlation tracking) which can lead to a 
sizable Doppler effect. At smaller angles, this 
motion is represented by the much smaller ra- 
dial velocity component ranging from 7 to 30 m/ s 
(|Williams fc PesnellllioTll iHathawav et al.ll2002t ). 
Our tests have shown that all main features of the 
GCI statistics discussed in this section, including 
the crossover scale £q, are clearly present in this 
central solar region where the supergranulation 
Doppler noise plays only a minor role. 

Unlike the autocorrelation (unipolar) dimen- 
sions, the a± dimensions of the quiet photosphere 
are systematically greater than 2 pointing at a 
consistently negative cross-correlation between op- 
posite polarity events. To visualize the cluster- 
ing geometries underlying the auto- and cross- 
correlation results, we constructed the unipolar 
and bipolar clustering maps defined respectively as 
M up = {vi-Y k \Pi = P k ] and M bp = {r, -r k \P t ^ 
Pk}- In essence, the maps are obtained by ad- 
justing the origin of the coordinate system conse- 
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quently with each of the N events, and superpos- 
ing the results in a single scatterplot while keeping 
track of the polarity of the events. As expected, 
the two maps are significantly different in shape 
(Fig. [5]). The dense core of the unipolar map 
represents the tendency of the magnetic events of 
the same sign to cluster into compact groups de- 
scribed by a± < 2. The empty central region of 
the bipolar map reflects negative correlations in- 
side the bipolar pairs leading to a± > 2, with the 
probability to find a neighbor of an opposite polar- 
ity increasing with the distance up to t ~ 30 — 40 
Mm. 



unipolar 




Fig. 5. — Unipolar and bipolar clustering maps of 
the quiet photosphcric region, M up and Mbp, vi- 
sualizing the difference in the geometry of the OE 
dynamics of the same (left) and opposite (right) 
magnetic polarity. 

The observed anti-correlation of the bipolar 
events in the quiet-Sun network implies a signif- 
icant imbalance of the apparent evolution of the 
magnetic flux in this region. Since the actual 
flux dynamics must be balanced (V • B = 0), 
our statistics manifest a geometric asymmetry of 



newly injected or canceling bipoles precluding the 
detection of a matching pole. It is known that 
an asymmetric coalescence of a preexisting un- 
resolved unipolar flux into magnetic concentra- 
tions that are sufficiently strong to be detected 
can produce unbalanced OEs: if two or more like- 
polarity ends of previously unresolved flux tubes 
come together, the average field strength can over- 
come the detection threshold for that polarity only 
( Lamb et al. 20081) . Likewise, asymmetric frag- 



mentation can lead to unbalanced demise events. 
Such sign-dependent effects are consistent with 
our analysis of the quiet solar region. We find that 
the positive events occurring in this region are by 
a factor of 1.5 - 3.0 (depending on r and B t h) more 
frequent than negative events indicating that the 
positive magnetic flux was significantly more frag- 
mented than the negative flux, which could natu- 
rally lead to the described artifacts. This interpre- 
tation is also supported by the fact that the dimen- 
sion a+ is systematically lower than a_ (see Table 
[T]) suggesting that spatial arrangement of positive 
polarity events along supergranule boundaries is 
more sparse compared to that of negative events. 

In addition to being a potential source of false 
event detection, the sign-dependent coalescence - 
fragmentation dynamics imply that on average, 
each negative event is magnetically connected with 
more than one positive event. The higher-order 
correlations associated with such magnetic struc- 
tures are beyond the scope of applicability of the 
proposed second-order statistical measures focus- 
ing on pairwise bipolar interactions such as those 
in magnetic loops. Addressing such correlations in 
future studies would be a challenging but rather 
important task. 

It is also noteworthy that in spite of the ge- 
ometric asymmetry of positive and negative po- 
larity fluxes, the origin and demise GCI dimen- 
sions of either magnetic polarity obtained for the 
quiet Sun are roughly the same. This can be a 
manifestation of a steady-state balance between 
the rate of injection of the magnetic flux from the 
convection zone, and the rate of its dissipation via 
coronal heating. It is generally accepted that the 
network magnetic fields play an important role in 
heating the solar corona. The total X-ray bright- 
ness of an active region is closely connected with 
the overall length of neutral lines and the pres- 
ence of strong magnetic shears in the neighboring 



13 



area (jFalconer et al.l I1997T) . Active core fields of 
the network of the quiet corona can be the main 
drivers of bulk coronal heating and dissipation as 
evidenced by the correlation between the location 
of bright X-ray e mission points and b ipolar mag- 



gested by our CGI measurements may play a vital 
role in initiating magnetic reconnection above so- 



netic structures ( Falconer et al.l 119981 ). Placed in 



this context, the nearly identical origin and demise 
GCI statistics strongly suggest that the coronal 
heating in the quiet Sun can be a magnetically 
driven process in which the sources and the sinks 
of the magnetic free energy share the same long- 
range correlation structure. 

The possibility of the evolving but statistically 
stationary quiet Sun network such as the one re- 
ve aled by our analysi s has been predicted earlier 
by I Simon et al. ( 2001 ) based on kinematic model- 
ing. Their results show that the total unsigned 
magnetic flux on the solar surface would decay 
within a few days due to the annihilation of oppo- 
sitely directed fields brought together by granular 
and supergranular flows, unless it is continually 
replenished at an appropriate constant rate. This 
requirement was presumably fulfilled in the stud- 
ied quiet solar region as suggested by the statisti- 
cally identical scaling behavior of the OE and DE 
dynamics. 

3.3. Active region NOAA 10365 

The correlation integral analysis of the active 
region leads to a rather different picture in which 
the origin and demise magnetic flux events follow 
distinct scaling laws, with a± « 2 for the former 
and a± < 2 for the latter events, see Fig. [4fc,d). 

The a± « 2 condition describing the OE popu- 
lation suggests the absence of sizable spatial cor- 
relations in newly detected bipolar pairs. This can 
be a signature of turbulent dynamics in the con- 
vection zone or some other process leading to a 
highly fragmented polarity separation line. Such 
lines are often associated with a sustained shear- 
ing motion of the field leading to a complex mix- 
ture of opposite magnetic polarities and resulting 
in the nonpotentiality of th e active region coronal 



field (|Schrijver et al.l 12005). Typically, the heat 



ing and flaring of coronal loops above an active 
region requires magnetic polarity inversion near 
at least one of the loop footpoints, and is greatly 
enhanced in the presence of a st rong magnetic 
shear in the core magnetic field (jFalconer et al.l 
19971 ). The disordered topology of the shear sug- 



lar active regions (lAntiqcho s 1998; P odgornv et al 



20071 iTorok et al.ll201ll) . 



In agreement with the proposed interpretation, 
strong shearing flows were observed ne ar the main 



polarity inversion line of NOAA 10365 (jChae et al 
I2004J ) . The region showed a systematic build- 
up of magnetic complexity en route to a flare- 
productive unstable configuration. The onset of 
coronal activity of NOAA 10365 was preceded by 
a pronounced distortio n of the shape of the main 
bipole separation line ( Verbeeck et al.ll201.il ). an 
injectio n of the helicity in the polarity inversion 
region ([Chae et aL 2004 ) , and a shar p increase of 
the Ising energy (I Ahmed et al. 2010) o f the bipo- 
lar magnetic flux ( Verbeeck et al.ll201ll ). A multi- 
wavelength analysis of NOAA 10365 confirms the 
formation of a highly twisted magnetic flux tube 
structure prior to M-class flares in the NOAA 
10365 region. Up to 50% of the helicity injec- 
tion rate was contributed by short-lived magnetic 
structures with ch aracteristic times 10-60 minutes 
(|Tian et al.1 120111 ) . The precursive process out- 



lined above could easily generate a rather complex 
photospheric field pattern with highly entangled 
magnetic polarities. In fact, the value a± ~ 2 im- 
plies that the complexity of the magnetic network 
of the active region was close to the theoretical 
limit in which the emergence locations of oppo- 
sitely directed flux elements are totally uncorre- 
lated. 

In contrast, the demise events in NOAA 10365 
are described by substantially lower a± values sig- 
naling positive correlations in bipolar magnetic 
pairs, possibly due to the cancellation of flux in re- 
connecting magnetic loops structures. This possi- 
bility is supported by the dependence of the cross- 
correlation dimension on the level of coronal flar- 
ing activity as shown in Fig. [6] The figure presents 
continuous scalograms a±(£,r) defined by eq. © 
characterizing the magnetic flux demise dynamics 
in flaring and non-flaring regimes of NOAA 10365 
(panels (b) and (c), correspondingly) as compared 
to the demise dynamics in the quiet photosphere 
(panel (a)). Each scalogram is accompanied by 
the r-dependence of small- and large-scale GCI di- 
mensions plotted in the same format as the right 
column in Fig. |4l 

The level of flaring activity is clearly encoded 
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Fig. 6. — Comparison of the correlation structure of the demise events the the quiet photosphere (a), in the 
NOAA 10365 active region without significant flaring activity (b) and in the same active region producing 
M-class flares (c). Left column: GCI exponents versus time scale r plotted in the same format as in Fig. 
ID Right column: continuous scalograms showing the dependence of the cross-correlation dimension a± on 
spatial and temporal separation between bipolar demise events. The red color prevailing in the flaring active 
region scalogram may be a signature of the magnetic flux cancellation. 
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in the shape of the scalograms shown in Fig. [6] 
The flaring state of the active region is character- 
ized by the condition a± < 2 (scalogram regions 
colored in red and yellow) for almost the entire 
studied range of i and r, confirming the presence 
of positively correlated bipolar demise dynamics 
at various spatial and temporal scales. The sign 
of the correlation changes from positive to slightly 
negative only in two narrow bands which could 
be associated with finite size effects imposed by 
the geometry of the active region. Even for very 
small delay times r ~ 1 min, there exist a distinct 
domain of I scales (~ 10 to 20 Mm) where the 
cross-correlation is positive. As was already men- 
tioned in Section 3.2, the quiet Sun exhibits no 
stable positive bipolar correlations. The narrow 
small-scale I band with a± < 2 seen in the top 
right panel of Fig. El can be an artifact reflecting 
insufficient number of closely located bipolar pairs 
in the quiet-Sun network. The quiescent active re- 
gion takes the intermediate position between the 
quiet Sun ans the flaring active region showing 
small but statistically significant positive correla- 
tions in selected I ranges. As in the case with the 
flaring NOAA 10365, the correlations describing 
nearly simultaneous demise events (r ~ 0) in the 
non-flaring state of this active region are only de- 
tectable short spatial distances. 

The observed tendencies suggest that the GCI 
dimensions of the active region indicative of pos- 
itively correlated bipolar demise events are con- 
nected with the flaring activity. In other words, 
there is a coupling between the symmetric reduc- 
tion of the photospheric flux and coronal dissipa- 
tion. The correlated demise events could be pro- 
duced by multiple bursts of magnetic reconnection 
in this NOAA 10365 leading to spatially localized 
symmetric reduction of oppositely signed line-of- 
sight magnetic fluxes. Indeed, the geometry of the 
chromospheric ribbons observed in the H a band 
unambiguously shows that the flaring activity of 
NOAA 10365 was triggered by magnetic reconnec- 
tion at coronal heig hts below the twisted flux tube 
of positive helicity ([Chandra et al. I l2009h . This is 
consistent with the results of 3D MHD simulations 
of the studied active region showing the formation 
of multiple current sheets in the vicinity of coro 
nal X-lines above NOAA 10365 (Podgornv et al 



mentary flare. The reconnection-driven annihila- 
tion of the coronal field above NOAA 10365 could 
lead to a detectable reduction of the photospheric 
flux in magnetically connected bipolar magnetic 
pairs, and contribute to the observed correlated 
occurrence of the demise events characterized by 
the condition a± < 2. To some extent, the bipo- 
lar GCI statistics can also represent simultaneous 
reduction of the magnetic flux in canceling bipolar 
pairs involved in bulk coronal heating around the 
active region (jSchriiver fc Titlell2002h . 

The t - dependence of a± dimensions of the ac- 
tive region is different than that of the quiet Sun, 
with positive correlations in disappearing bipolar 
pairs being the strongest at smallest r scales. This 
suggests that the characteristic time scale of the 
apparent flux cancellation can be comparable with 
the 1-minute image sampling time. The correla- 
tions also increase with the decrease of the spatial 
scale I and tend to be less pronounced at I > £ - 
If the association of these events with canceling 
magnetic flux elements is correct, our observations 
indicate the dominant contribution of small- to 
medium size coronal loops in the reconnection pro- 
cess. 

Positively correlated unipolar demise events 
are reliably identified for B t h exceeding ^2 stan- 
dard deviations of the active region field but are 
not seen at lower detection thresholds. There- 
fore, they are a property of an intense magnetic 
field. For higher thresholds, the cross-correlations 
monotonously increase with B t h suggesting an in- 
volvement of non-Gaussian magnetic field fluctu- 
ations, such as those associated with intermittent 
turbu lence below the photosphere ( Abramenkol 
20051 ) or transient multisc ale dissipation in the 



corona 



Uritskv et al] ([20071 ). This explanation is 



in line with earlier observations showing that the 
preferred locations of bright coronal points corre- 
late wit h the regions of intense unsigned magnetic 



flux (see lBenevolenskaval (|2010h and refs therein), 



2007; IPodeornv fc Podgornvl 120081 ) . Each of the 



current sheets may well be responsible for an ele- 



although the physical mechanism of this depen- 
dence likely to be driven by topological transitions 
in a complex multiscale magnetic network remains 
to be understood. 

4. Conclusions 

Our findings demonstrate that the developed 
data analysis approach based on the GCI algo- 



16 



rithms and direct spatiotemporal tracking of mag- 
netic elements is capable of addressing various as- 
pects of the dynamic of solar photospheric net- 
work, including the geometry of supergranular 
structures, emergence of the nonpotential flux, 
and cancellation dynamics of bipolar magnetic 
pairs. 

The results obtained reveal a significantly dif- 
ferent multiscale dynamics of the photospheric 
magnetic flux in quiet and active solar regions 
as represented by the studied sample of SOHO 
MDI images. Quiet Sun exhibits robust unipo- 
lar clustering of both OE and DE events over a 
broad range of photospheric distances correspond- 
ing to meso- and supergranulation, with a small- 
scale quasi one-dimensional clustering along su- 
pergranular boundaries and a more isotropic clus- 
tering at the level of supergranules. These clus- 
tering regimes place important constraints on the 
topology of the photospheric network and spa- 
tially distributed storage of free magnetic energy 
in the quiet Sun as illustrated by Fig. [TJa). No 
detectable cross-correlations in bipolar magnetic 
pairs which could be related with unstable mag- 
netic loops were found in the quiet photospheric 
region. NOAA 10365 solar active region exhibits 
signatures of non-potential magnetic structures 
with a topologically complex polarity separation 
line. 

Flaring activity in the active region is accom- 
panied by positively correlated magnetic demise 
events which can be a signature of magnetic re- 
connection. A possible underlying geometry which 
would lead to the observed CGI signatures of the 
active region is a sheared loop arcade with an ir- 
regular polarity inversion line as depicted in Fig. 
[7jb,c). The topological complexity of this struc- 
ture can be instrumental because simple shearing 
of an arcade is typically insufficient to any signif- 
icant manifestati on of solar actiyity su ch as flares 



or prominences (jDahlburg et al.1 119911 ). The sta 



tistical results obtained for AR 103675 speak in 
favor of so-called sympathetic eruptions in this 
complex active region. Such eruptions may oc- 
cur within a relatively short period of time in 
spatially separated flux ropes due to a variety of 
MHD-scale trigger s destabilizing the coronal field 



(jTorok et alj |2011) 




Fig. 7. — Examples of spatial arrangement of mag- 
netic events in the quiet Sun (a) and in the active 
solar region in quiescent (b) and flaring (c) states 
consistent with our analysis results. Red lines in- 
dicate positively correlated demise events of oppo- 
site magnetic polarity presumably associated with 
magnetic flux cancellation events in the active re- 
gion whic h could be involved in sympathetic solar 
eruptions (jTorok et al.ll201ll) . 



tical problems associated with solar magnetic 
complexity and its impact on solar activity. For 
example, is there a relationship between statis- 
tical properties of fractality/complexity in the 
photosphere and corona? Some studi e s pro - 



vide definitive answe r (e.g.. IConlon et al.l (120101): 



Falconer et al. ( 20091 ); Abramenko fc Yurchvshvn 
( 2010a H)), but some others argues ag ainst such 



The proposed methodology can be helpful in 
resolving a broader set of theoretical and prac- 



possibility (|Dimitropoulou et al.l 120091) . If the 
complexity linking between the two solar struc- 
tures exists, one can expect the CGI signatures 
of photosphere and corona to vary in a consistent 
way. 

It is also not clear how the fundamental plasma 
processes responsible for the formation and evo- 
lution of the magnetic network organize them- 
selves across such a broad range of scales. Ac- 
curate evaluation of CGI dimensions for statis- 
tically representative sets of solar regions and 
conditions may provide important clues on the 
physics of the underlying process, for instance, 
by validating certain typ es of turbulent cascades 
ijMuller fc Biskamdl2000h . The proposed method- 
ology can also help to better understand how the 
scaling properties of the quiet and active Sun net- 
works relate to each other, and how these prop- 
erties reflect the formation of unstable magnetic 
topologies leading to plasma heating and dissipa- 
tion. The relationship between the temporal and 
spatial scaling laws governing multiscale geometry 
and dynamics of solar magnetic field is another im- 
portant subject of future research. Scaling charac- 
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teristics of sheared magnetic structures generated 
by oppositely directed line-of-sight flux elements, 
and the contribution of these structures the coro- 
nal dissipation in quiet and active solar regions are 
also largely unexplored and requires careful con- 
sideration. To clarify the underlying physics of 
the observed effects, it would be essential to com- 
bine our statistical tools with conjugate observa- 
tions o f coronal act ivity as well as solar tomog- 
raphy (|Davilalll994l) and other three-dimensional 
reconstruction techniques. 

The developed approach can also be used 
for improving the forecasting of flaring activity 
through identification of topologically complex 
unstable magnetic configurations. Spatial mag- 
netic complexity of an active region as measured 
by the power-law decay of the Fourier spectrum of 
the longitudinal field was shown to c orrelate with 
the fl are productivity of the region ( Abramenkol 
20051 ). This dependence suggests that active re- 
gions with more intermittent magnetic structure 
tend to produce higher energy flaring events. Our 
methodology offers a more comprehensive way of 
measuring magnetic intermittency characterizing 
each magnetic polarity (based on a + and a_), as 
well as bipolar magnetic structures (a± ) . Analysis 
of these correlation measures, as well as the tem- 
poral CGI dimensions not considered here, may 
help better understand photospheric drivers of 
bulk coronal heating and local topological precur- 
sors of bright coronal points and flaring activity 
across the entire range of the involved spatial and 
temporal scales. 
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